Emergence of canonical ensembles from pure quantum states 
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We consider a system weakly interacting with a bath as a thermodynamic setting to establish a 
quantum foundation of statistical physics. It is shown that even if the composite system is initially in 
an arbitrary nonequilibrium pure quantum state, the unitary dynamics of a generic weak interaction 
almost always drives the subsystem into the canonical ensemble, in the usual sense of typicality. 
A crucial step is taken by assuming that the matrix elements of the interaction Hamiltonian have 
random phases, while their amplitudes are left unrestricted. 
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The concept of typicality has recently gained renewed 
interest in the context of the quantum foundation of sta- 
tistical physics [H-Q. In particular, the following has 
been shown [l|: Consider a closed quantum system hav- 
ing a very large Hilbert space dimension. One randomly 
picks a generic pure quantum state among those having 
a particular energy within an energy shell. It then turns 
out that the state of a small subsystem with the rest be- 
ing traced out almost always, i.e., typically, has a canon- 
ical distribution. Conceptually, this perspective seems to 
contrast with the traditional one, in which the canonical 
ensemble of a subsystem follows when the total system 
is in the maximum-entropy state, i.e., the microcanoni- 
cal ensemble jij. In the new perspective, the entropy of 
the total system seems to be zero, while its subsystem 
nevertheless has a large entropy. 

It should be noted here that the above two seem- 
ingly different perspectives are equivalent in the following 
sense: The new perspective implicitly assumes ergodic- 
ity as the pure state is taken with respect to the uniform 
measure. The initial set of states from which one is taken 
is thus, due to symmetry, a microcanonical ensemble. 
As the objective realism is an inherent nature of those 
statistical descriptions and the microcanonical ensemble 
is unitarily invariant, whichever basis is chosen, picking 
one state out of the ensemble, having prior and complete 
knowledge of the system and hence making the entropy 
vanish in advance, should not make any difference to the 
physical consequences. An important point is the fact 
that when we are given a thermodynamic system, we are 
already observing only one state among those in the en- 
semble. Whether we know beforehand which state it was 
does not matter. What statistical physics asserts is that 
we indeed need not know such microscopic information 
because almost all individual states will yield the same 
statistical quantities as long as their fluctuations vanish 
in the thermodynamic limit. 

It is thus instructive to discriminate between two as- 
pects behind the typicality arguments [lj. On one hand, 
the typicality is generally a property of a large Hilbert 
space dimension. On the other hand, its actual manifes- 
tation (e.g., the resulting subsystem state) depends on 



the choice of the initial ensemble from which the state 
is taken. Hence, there remains a fundamental question 
yet to be answered: Why do we pick a state with equal 
priority? Why does the system have no preference to par- 
ticular states? Recall that the ergodic hypothesis is not a 
necessary physical consequence that every system should 
result in, hence, hypothesis. A more relevant question 
should be when and how the microcanonical ensemble is a 
plausible representation of a given system, which we shall 
consider in this paper. Whereas classical concepts such 
as chaos and mixing support the ergodic hypothesis on 
a mathematically rigorous footing, there are no stringent 
counterparts of them in quantum physics, which poses 
difficulties Note, however, that ergodicity and mixing 
are empirically defined, not derived from first principles, 
and by nature statistical concepts defined on a measure 
space. In order to establish a corresponding quantum 
theory, the question is, as for any statistical argument, 
how to introduce some statistical nature, i.e., random- 
ness, in a physically plausible manner. In this regard, it 
is tempting to consider the inherent phase ambiguities in 
quantum systems as a natural source of randomness. 

In this paper, we consider a closed thermodynamic sys- 
tem, composed of a system S and a bath B, in an arbi- 
trary nonequilibrium pure quantum state with a certain 
energy, and study how it evolves when a generic weak in- 
teraction between S and B, whose meaning is clarified be- 
low in terms of phase randomness, is brought in. We show 
that for almost all interactions system S is driven to the 
canonical ensemble under some reasonable conditions. 
To this end, we separate the time-independent Hamilto- 
nian H into the noninteraction part Hq — H$ + Hb with 
the system (bath) Hamiltonian Hs (Hb) and the inter- 
action part Hi. We expand the state at time t with the 
eigenstates \j) of H as \^{t)) — J2j c j(t) The tran- 
sitions between the energy levels are determined by the 
off-diagonal elements Hju = (j \ Hi \k) (we let Hjj = 0). 

Our starting point is to note that, given the inherent 
complexity and arbitrariness of thermodynamic systems, 
it is practically reasonable to assume that the phases of 
Hjk, which vary sensitively with boundary conditions, 
positions of interaction, etc., are determined randomly. 
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Let us denote by Ojk the phase of Hjk- Our basic as- 
sumption is that Ojk £ [0, 2ir) could be taken indepen- 
dently and randomly with respect to the uniform mea- 
sure. For later convenience, however, we transform the 
set of 6jk into an equivalent set of uniform random vari- 
ables <j>j,0'- k £ [0, 2n) for j, k ^ 1, where cf)j — Oji and 
9'- k = Ojk — {(f) j — 4>k) (let us set (f>j = 0'- k = otherwise). 
We then take only <f)j as the random variable while letting 
e l9 ik De absorbed into Hjk, which is enough to account 
for the typicality of thermodynamic evolutions. We now 
have 

H jk -> jfo-MH jk , 4>j,^k G [0,2tt). (1) 

In what follows, we shall denote by (■} the average over 
the ensemble defined with respect to the random vari- 
ables 4>j Q. This ensemble reminds us of the random 
matrix theory Q. Our assumption is, however, much 
weaker because except for the phase randomness, no as- 
sumptions are made on the statistical distribution of their 
magnitudes. The formalism below also holds, with slight 
modification, for real Hjk with random signs. 

As discussed before, a source of randomness is a nec- 
essary element of any statistical argument. In this sense, 
we simply take the hypothesis of Eq. (fTJ) for granted with- 
out a supporting microscopic theory, except for argu- 
ing that there is no physical law for an arbitrary sys- 
tem to prefer a particular phase unless specially pre- 
pared. We are just concerned with its consequences 
when such an ensemble is acceptable as properly describ- 
ing a given system. Let us take the interaction picture 
to describe the dynamics in terms of the time evolu- 
tion operator U(t,t ) = 1 -i f* Vi(t')U(t',t )dt' with 
Vi(t) = e iHot H ie - iHot . The coefficients Cj(t) can be 
written as 

to)ck{to) (2) 

with Uj k {t,t ) = (j\ U(t,t ) \k). From Eq. ([1]), we have 
(U* k (t,t )U lm (t,t )) = \U jk (t,t Q )\ 2 S 3l 8 km (3) 

with the Kronecker deltas 8j k and 8km- Note that the 
right-hand side is invariant under the phase randomness. 
Eq. ([3]) implies that 

(c*(t)c k (t)) = (\c 3 (t)\ 2 )S 3k (4) 

is satisfied at every t. Note that this condition along with 

(\ c j(t)\ 2 ) — constant (5) 

constitutes the defining character of the ergodic state, or 
equivalently the microcanonical ensemble. Once condi- 
tion Q is understood, the remaining question is how to 
derive condition ([5]) from it. 



Before proceeding, we clarify three extra (but rea- 
sonable) conditions. The first condition is concerning 
the time scale. Among various time scales characteriz- 
ing dynamical systems Q, two time scales are relevant 
in our context. The first one is the Heisenberg time 
t~h ~ 1/ AE, where AE is the characteristic energy level 
spacing of the system. This time scale is related to the 
energy-time uncertainty. When t <^_th, the system can- 
not see the discrete nature of its energy spectrum. In 
this case, if the dynamics is diffusive, it persists as in 
the classical counterpart system. On the other hand, 
when t > th, the system recognizes the discreteness and 
the dynamics can deviate from the classical counterpart. 
Dynamical localization Q , a dynamical analog of Ander- 
son localization , is a well-known phenomenon in this 
regime, demonstrating a stark difference between quan- 
tum and classical dynamics. Another relevant time scale 
is the ergodic time te, which estimates the time for the 
system to wander the entire region of the available phase 
space. This time scale is largely dependent on the mi- 
croscopic details, and discussing it quantitatively is out 
of our scope. Qualitatively, we can argue that the more 
complex the transition structure is, the faster is the dif- 
fusion process, hence the shorter is te- Later it becomes 
evident that one can introduce such a concept as con- 
nectivity in the purely classical sense. Here we simply 
make an inevitable assumption that the transition struc- 
ture is such that all the energetically allowed levels are 
indeed connected through finite transitions (in the classi- 
cal sense) and its time scale meets te <C t#. If te > th, 
dynamical localization can occur, inhibiting the system 
from reaching the ergodic state. Given the fact that the 
total energy increases linearly with the number of par- 
ticles, while the Hilbert space dimension increases expo- 
nentially, for ordinary thermodynamic systems, the en- 
ergy spectrum is continuous in any physically realistic 
sense and the Heisenberg time can never be reached. The 
condition te *C th is thus automatically met in all such 
systems as long as the transition structure is as assumed 
above. At this point, it is worthwhile to compare ear- 
lier dynamical considerations in the present context with 
ours [2j. Therein, the reason for being able to observe 
an ensemble behavior was ascribed to the fact that the 
time average of the relative phase (e^ l ^ E ^ Ek ^ between 
two energy eigenstates effectively vanishes in a long time. 
However, this time averaging effect is relevant after the 
Heisenberg time. In the thermodynamic systems we con- 
sider, this effect cannot be observed. 

The second condition is the high dimensionality, which 
is essential in statistical physics. We quantify this in 
terms of the effective dimension, or the effective num- 
ber of involved states, do(t) -1 = Ylj l c i( f )| 4 ; which es- 
timates how many coefficients are nonzero. We assume 
da(t) 3> ds is met at every i, where d$ is the Hilbert 
space dimension of system S. As our regime is where the 
total system has an extremely denser energy spectrum 
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than system S, this condition is also automatically met. 

The last condition is concerning the interaction be- 
tween the system and the bath. As the spectrum is 
a quasicontinuum, we can write the state as \j) — > 
\Es,Eb), where Es and Eb denote the energies of 
system S and bath B, respectively. We assume that 
\{Es,Eb\H\\E' s ,E' b )\ is a quasicontinuous function of 
these energies, or at least of Eb and E' B (recall that the 
bath has a much denser energy spectrum). For example, 
this is the case in scattering problems, which would de- 
scribe the dynamics of gases [llj. In addition, we assume 
that the interaction is weak. This has a clear analog in 
statistical physics, wherein the derivation of the canoni- 
cal ensemble from the microcanonical ensemble relies on 
the assumption of weak interaction. 

We are now in a position to proceed to our main result. 
From Eq. @, the population in \ j) can be written as: 

pm = \cM\ 2 =p%t) + ip f /(t) +p } r(t) +p f /(t)}, (6) 

where 

P %t) Ht/„(Mo)l 2 MMI 2 + ]T \u 3k {t,t )\ 2 \c k (t )\ 2 , 



Pi 2 (*)= E MMo)tf*KMo)cfc(ioK(to)(l-<5, 



(7) 
(8) 



(9) 



Our plan is to calculate the average dynamics, wherein 
system S asymptotically reaches the canonical ensemble 
by tracing out over bath B. We then show its fluctuation 
vanishes, which proves that for almost every individual 
dynamics system S reaches the canonical ensemble. 
From Eq. it is easily seen that 



(j> j (t)) = (p<*(t)). 
By defining the transition rate 

Wk^j(t,t ) = \U jk (t,t )\ 2 



(10) 



(11) 



and using f t \Ujj(t, t )\ 2 = - f t \U kj (t,t a )\ 2 , which 

follows from the unitarity of U(t,to), Eq. ([TO)) reads 



dt 



+ Y2 w k^j( t , t o) (Pk(*Cl)) 
kjtj 



(12) 



As the interaction is weak and the energy spectrum is 
continuous, the transition rate is given by the Fermi 
golden rule: 

w k ^ 3 = 27r|(j| H t \k)\ 2 8{Ej - E k ), (13) 



where it is understood that the density of energy levels 
T(Ej) appears when the summation is taken over the fi- 
nal states around \j), and Ej and E k are the correspond- 
ing energies [l2j]. Eq. (fT3|) can also be calculated up to 
higher orders using the so-called T matrix, but the result 
is essentially the same in nature 
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For a continuous energy spectrum, transition should be 
treated carefully with respect to a well-defined normal- 
ization method (l2j |. In our case, the Hilbert space can 
be naturally separated into ds subspaces, H = 0„li H^, 
where is defined by a projector (^^| ® I B with 
|^^) the juth eigenstate of H$ and I B the identity op- 
erator for bath B. We are concerned with the transitions 
between these subspaces, which represent the transitions 
in system S of our interest. 

Suppose the total system has energy E within a cer- 
tain precision, including that due to the energy-time un- 
certainty. Let us define the population in "H M : 



(14) 



\j)en„ 



We assume Wj^ k = for |j) , |fc) S H^, since such transi- 
tion within a subspace can be included in Hb ■ From the 
definition of H M , it can be seen that the density of energy 
levels of is solely determined by that of the bath. De- 
noting by E^ the eigenvalue of j^f), it can be written 
as T B (E - E$ ) . Note that the transition of Eq. (fTOj) pre- 
serves energy. By suitably summing Eq. (|12l) and taking 
into account the energy spectrum being continuous, we 
can obtain 



UjtfJ, 



E 



,T B {E-E s u )(P^t )) 



l^b(E — E^) (P u (to)) , 



(15) 



where W^ v = 



denotes 



Wa 



= 2ir\(E^E- 



E^H^E^E 



(16) 



It is instructive to represent this situation as a transition 
network of ds nodes. In particular, as the transition rate 
is independent of time, the transition structure is visu- 
alized as a Markov chain. One can define a population 
vector P(t) from (P^(i)) and a transition matrix T such 
that 



J2 W ^ r B(E-E B )St, 
W^T B (E-E B )St. 



The population is then given by 

P(t + n6t) =T n P(i ) 



(17) 
(18) 

(19) 



for a time interval St chosen to be fairly small, but large 
enough for the Fermi golden rule to be valid. 
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This picture is very useful because the ergodic theory 
is well established in the Markov chain . When every 
node is reachable from everywhere, the Markov chain is 
called ergodic or irreducible. In this case, the ergodicity 
means in the physical context that the time average of the 
population approaches a steady state in that it is given by 
linin-^oo E"=o P(jdt)]/( n +l)- An ergodic Markov chain 
is called regular if some power of the transition matrix 
has only positive elements. A regular Markov chain has 
a stronger meaning of ergodicity in that P(t) itself ap- 
proaches the steady state as t — !> oo. In either case, there 
is a unique steady state determined by detailed balance. 
From Eq. (|15l) . the steady state is found to be 

(P„(t)><xr fl (20) 

Note that this condition is equivalent to that of the 
canonical ensemble, where it is understood Tb(E) in- 
creases exponentially with E Q • A similar result can be 
obtained using the master equation approach |14j . 

Finally we show that the fluctuation of the average 
dynamics effectively vanishes. Note that Pj(t) in Eq. ([7]) 
is invariant under the phase randomness for any given 
state at time to- The off-diagonal contributions Pj(t) 

and pj (t) in Eq. ([6]) can thus be interpreted as random 
fluctuations to Pj(t). It is thus enough to show that the 

variances of pf 1 (t) and p^ (t) vanish in our regime. These 
are calculated straightforwardly from Eq. ^ . We have 

(|pJ 1 W| 2 > = (l^'(t^o)| 2 |c J (i )| 2 El f/ ^(*'*o)| 2 |c fc (io)| 2 ) 

<yjlMMo)| 2 (Mio)| 2 )- (21) 

The last line follows from the Cauchy- 
Schwarz inequality (J2k Wjk{t, *o) | 2 |cfc(*o) | 2 ) < 
(Efcl^'fc(*>*o)| )(EfcM*o)| ) and the property of 
unitary operators Ylk Wjk(t> to)\ < 1- hi the same 
fashion, the variance of pj (t) is calculated as 

(\p f ;(t)\ 2 ) < /r£|[/; fe (Mo)| 2 (Mto)| 2 ). (22) 

By summing Eqs. (|2ip and (|22p. we have 

<|AP M (t)| 2 ><2y^(P M (*)). (23) 

For any statistical distribution {P^(t)} of ds variables, 
this variance is negligible since do is exponentially larger 
than ds- 

In summary, we have presented a quantum descrip- 
tion of the generic thermalization process for a subsystem 
weakly interacting with a bath. A crucial step was taken 



by assuming that the matrix elements of the interaction 
Hamiltonian have random phases as in Eq. ([1}, while 
their amplitudes are left unrestricted. Granting three 
extra conditions stated after Eq. ([5]), we found that the 
average subsystem dynamics reduces to a Markov chain 
of Eq. ([H?f whose steady state coincides with a canonical 
ensemble of Eq. (|2"0)l . Indeed we found that such an av- 
erage dynamics is the actual manifestation of almost ev- 
ery individual dynamics as its fluctuation vanishes. This 
implies that for an arbitrary initial state and a generic 
unitary dynamics of the whole system as characterized 
above, the subsystem is almost always led to a canonical 
ensemble. The irreversibility is not problematic because 
our time scale is such that the recurrence in quantum 
dynamics is not observable [l^. 
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